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ABSTRACT 

The quasilinear theory for cosmic ray propagation is a well known and widely accepted 
theory. In this paper, we discuss the different contributions to the pitch-angle Fokker- 
Planck coefficient from large and small scales for slab geometry using the damping 
model of dynamical turbulence. These examinations will give us a hint on the limita- 
tion range where quasilinear approximation is a good approximation. 
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1 INTRODUCTION 

The propagation of cosmic rays (CRs) is affected by their 
interaction with a magnetic field. This field is turbulent 
and therefore, the resonant interaction of cosmic rays with 
MHD turbulence has been discussed by many authors as the 
principal mechanism to scatter and isotropize cosmic rays 
(Schlickeiser 2002, Yan & Lazarian 2002). Although cosmic 
ray diffusion can happen while cosmic rays follow wander- 
ing magnetic fields ( Jokipii 1966) , the acceleration of cosmic 
rays requires efficient scattering. For instance, scattering of 
cosmic rays back into the shock is a vital component of the 
first order Fermi acceleration (see Longair 1997). 

One important interaction between CRs and MHD turbu- 
lence is gyroresonance scattering. While the propagation of 
most moderate energy CRs don't sense slower time varia- 
tions of turbulence, turbulence is dynamical for low energy 
CRs or when the parallel speed of CRs becomes comparable 
with Alfven speed v« ~ Va- Dynamical MHD turbulence 
leads to resonance broadening so that a Breit-Wigner-type 
function instead of a S function should be used (Schlickeiser 
& Achatz 1993, Bieber et al. 1994). In this case CRs interact 
with a whole range of turbulence from large to small scales. 

Scattering of CRs cannot happen by the magnetic wave 
which frequency in the particle frame is much lower than 
the Larmor frequency. Indeed if magnetic field changes much 
slower than the particle gyrorates, the quantity v\/B is pre- 
served (Landau & Lifshitz 1957). As the result while the 
pitch-angle of the particle changes as the wave passes by 
the total change of the angle is zero. 

In quasilinear theory (QLT) the assumption of unper- 
turbed orbit results in non-conservation of the adiabatic in- 
variant p = mv\/2Bo. Whereas small scale contribution 
corresponds to the sharp resonance in magnetostatic limit, 



the contribution from slow large scale can be overestimated 
by QLT since the adiabatic invariant ought to be conserved 
when the electromagnetic field varies on a time scale longer 
than the gyroperiods of CRs (Chandran 2000, Yan & 
Lazarian 2003). 

In QLT the turbulent field is presented in terms of Fourier 
modes. It is clear that the slow (compared to the particle 
Larmor frequency) components cannot scatter cosmic rays 
as the consequence of the preservation of the adiabatic in- 
variant. However, this effect is not treated in the conven- 
tional QLT, which has been widely used to obtained a lot 
of results mostly within the framework of the slab model of 
turbulence. How reliable are those results? We attempt to 
answer this question within this paper. 

Therefore, we will discuss the different contributions 
from large and small scales for the slab model of turbulence. 
We will give the regimes where the interaction is dominated 
by large and small scale contributions. This will give us a 
hint on the limitation range where quasilinear approxima- 
tion is a good approximation. Because QLT is considered as 
a standard tool for the calculation of diffusion coefficients 
(e. g. Bieber et al. 1994, Chandran 2000, Schlickeiser 2002, 
Yan & Lazarian 2003) it is of great significance to explore 
the validity of the quasilinear approximation. 

Within quasilinear theory the parallel mean free path Ay re- 
sults from the pitch-angle-cosine (fi = p\\/p) average of the 
inverse of the pitch-angle Fokker-Planck coefficient as 
(Jokipii 1966, Hasselmann & Wibberenz 1968, Earl 1974) 
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The pitch-angle Fokker-Planck coefficient is calculated from 
the ensemble-averaged first-order corrections to the particle 
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orbits in the weakly turbulent magnetic field (Hall & Stur- 
rock 1968) 



ZV(^) = Re I di< (i(t)fi*(t + 0> 



(2) 



and depends on the nature and statistical properties of 
the electromagnetic turbulence and the turbulence-carrying 
background medium. 

In Sect. 2 we derive general expressions for the pitch-angle 
Fokker-Planck coefficient D MM . In Sect. 3 we compare the 
small scale and large scale contribution to the Fokker-Planck 
coefficient with each other. This results can be used to cal- 
culate the cosmic ray parallel mean free path using Eq. Q 
what is done in Sect. 4. 



2 THE PITCH-ANGLE FOKKER-PLANCK 
COEFFICIENT £> MM 

The key input to the calculation of Fokker-Planck coeffi- 
cients is the correlation tensor Py. It is common to assume 
the same temporal behaviour for all tensor components 

p i3 -(M) = Py(£,o)-r(M) (3) 

where F is the dynamical correlation function. In the past 
several models for F were discussed (Schlickeiser 2002). In 
the current paper we use the damping model of dynamical 
turbulence (DT-model, Bieber et al. 1994) for slab geometry: 



r(fc,t) = r(fc||,t) = e 
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Here va is the Alfven speed and a is a parameter which 
allows us to adjust the strength of dynamical effects, a = 
corresponds to the magnetostatic limit whereas a = 1 
corresponds to strongly dynamical turbulence. 

A further important input to our calculations is the 
wave spectrum. For the following examinations we use a sim- 
ple power spectrum without energy- and dissipation-range: 
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Here and in the following discussions SB is the total strength 
of the turbulent magnetic field and Bo is the strength of 
the magnetic background field (mean field). The parame- 
ter s in Eq. (JSJ is the inertial-range spectral index. For a 
Kolmogorov-spectrum we have s = 5/3 which was consid- 
ered by Teufel & Schlickeiser 2002 and 2003. Another choice 
would be s = 3/2 (see Cho & Lazarian 2002). Assuming 
this power spectrum, slab geometry and the DT-model the 
pitch-angle Fokker-Planck coefficient can be written as (see 
Teufel & Schlickeiser 2002) 
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Now we split this integral at the wavenumber k = ek r 

(s-l)Q. 2 (l- n 2 ) (8B\ 2 
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The resonant wavenumber k re3 is defined through (see Yan 
& Lazarian 2002) 
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where Rl is the gyroradius. In the first integral we have 
fc _1 e _1 /iJ?L and the inverse wavenumber is larger than 
several gyroradii if e <C 1. We call this case the large scales. 
In the second integral we have fc _1 ^ e _1 /iRl and the invers 
wavenumber is smaller than several gyroradii. We call this 
case the small scales. We always have 
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where we used the dimensionless rigidity R = Rhk m in- In 
the current paper we restrict our analysis to medium rigidi- 
ties {R <C e <C 1, v ^> ova) and we have therefore 



< e ■ k r 



(11) 



In this and the following equations we use fj, =\ fj, |, ^ a ^ 
1 and 1 < s < 2. With 
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and the small scale Fokker-Planck coefficient 
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Now we split the large scale Fokker-Planck coefficient into 
two parts 
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and use the transformation x — k m i n /k in the first and 
x = ek rf:s /k in the second integral to obtain 
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(17) 
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where we defined the both integrals: 
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with 

a= — >1. (19) 

We can also use the transformation x — ek res /k to simplify 
the small scale Fokker-Planck coefficient 



The total Fokker-Planck coefficient is equal to 
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2.1 Analytical results 

As demonstrated in Teufel & Schlickeiser 2002 the both in- 
tegrals list can be calculated approximatelly for different 
cases. For medium rigidities we have \iR <C e <C 1 and 
a/J?> 1 and we find: 
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With these approximations we obtain the following expres- 
sions for the total, small scale and large scale Fokker-Planck 
coefficient: 
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The magnetostatic result (a = 0) for the total pitch-angle 
Fokker-Planck coefficient is 
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for all values of /i. 



2.2 Numerical results 

To test the analytical predictions we also calculated the dif- 
ferent Fokker-Planck coefficients numerically. To calculate 
D MM or A || as a function of the dimensionless rigidity R we 
must express the ratio va/v through the parameter R: 
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Fig. and [5] show the different Fokker-Planck coefficients 
for protons and the following set of parameters: 
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(27) 

With these parameters we have Ro(p + ) ~ 0.169 and 
Ro(e~) ~ 92.1 ■ 10 -6 . Fig.^shows the numerical calculated 
Fokker-Planck coefficients for e = 1. In this case we split 
the integrals at k = k res . As expected the small scale and 
and large scale Fokker-Planck coefficients are approximatelly 
equal except for small a where the large scale contribution 
is dominant. In Fig. [5] we have shown the numerical results 
for e = 1/3. Here the small scale contribution is always dom- 
inant except for small /i. 



3 COMPARISON OF THE SMALL SCALE 
AND LARGE SCALE FOKKER-PLANCK 
COEFFICIENTS 

In this section we calculate the ratio 
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for the cases a/i -C e -C 1, e <C a/t <C 1 and e < 1 < a/i. 

3.1 The case a/i < e < 1 
Here we find with Eq. @ 
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Figure 1. Numerical results for the different Fokker-Planck coef- 
ficients for protons for e = 1. The solid line shows the total Fokker- 
Planck coefficient in comparison with the magnetostatic results 
(dotted line), the small scale Fokker-Planck coefficient (dashed 
line) and the large scale Fokker-Planck coefficient (dash-dotted 
line). 



3.2 The case e < ap. < 1 

In this case we have 
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and therefore 
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For e <C an <C 1 the large scale contribution is much smaller 
then the small scale contribution. 



3.3 The case e < 1 < afx 

Here we find with Eq. (23 or Eq. (23 
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Also in the case e < 1 < a/i the large scale contribution is 
much smaller then the small scale contribution. 
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Figure 2. Numerical results for the different Fokker-Planck co- 
efficients for protons for e = 1/3. The solid line shows the to- 
tal Fokker-Planck coefficient in comparison with the magneto- 
static results (dotted line), the small scale Fokker-Planck coeffi- 
cient (dashed line) and the large scale Fokker-Planck coefficient 
(dash-dotted line). 



£{aa < e < 1) = 

and therefore 
D LS 3> D ss 

MM " MM " 



2sin(7rs/2) 



— | »1 

au 



(29) 



(30) 



The last both results can also be derived directly from Eq. 
12311 . In the case a/i C e < 1 the large scale contribution is 
much larger then the small scale contribution. 



4 THE PARALLEL MEAN FREE PATH 

In this section we calculate the parallel mean free path us- 
ing the results for the Fokker-Planck coefficients of the last 
section. 



4.1 The total contribution 

With Eq. (23 we have 
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which is approximately equal to the magnetostatic results 
(A|^ s ). For a Kolmogorov spectrum (s = 5/3) we find that 
the parallel mean free path is ~ R 1 ? 3 , whereas for s = 3/2 
the parallel mean free path is ~ R 1 ' 2 . 

We also find that the result for Ajj°' is independent of the 
damping parameter a. Therefore we come to the conclusion 
that dynamical effects can be neglected if we calculate the 
total parallel mean free path for medium rigidities and for 
slab geometry. This result is in agreement with numerical 
calculations presented in this paper (see Fig. and [3 ■ 



Spurious contribution to CR scattering calculations 5 



4.2 The small scale contribution 

If we neglect the large scale contribution if we calculate the 
parallel mean free path we have 

ss _ 3, f 1 (1-m 2 ) 2 



3v 
4 

3d 



l/a 







r Dfiiian < e < 1) 



(1"M 2 ) 2 



+ — / d M 

' l/a 

3w Z" 1 (l~/x 2 ) 2 
+ 4 J 1/e dfX D™(e«l«a„) 



(37) 



and together with Eq. 1231 we find for the parallel mean free 
path 
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If we neglect the large scale contribution the parallel mean 
free path goes to infinity. This result comes from very small 
values of \i. It should be noted that in this regime QLT is no 
longer valid because for fi — > nonlinear effects are essential 
(see Shalchi 2005). If we would use a more accurate, non- 
linear theory instead of QLT, the problem of the infinitely 
large parallel mean free path (Eq. 1381 ) should vanish. 



4.3 The large scale contribution 

If we neglect the small scale contribution if we calculate the 
parallel mean free path we have 
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With Eq. >ti we obtain for the parallel mean free path 
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In this case the parallel mean free path is larger than the 
magnetostatic result. This would be the case when the tur- 
bulence is damped at a scale larger than the resonant scale. 



5 DISCUSSIONS AND SUMMARY 

In the current paper we calculated the pitch-angle Fokker- 
Planck coefficient and the parallel mean free path using 
QLT. Our intension was to find out how important dynami- 
cal effects are for pure slab geometry and for medium rigidi- 
ties and whether the small scale or the large scale contri- 
bution controlls the total Fokker-Planck coefficient and the 
parallel mean free path. 

If we consider Eq. I36H we come to the conclusion that 
we can neclect dynamical effects if we calculate the parallel 
mean free path for medium rigidities. The situation would 



be different if we would consider small rigidities and if we 
would use a power spectrum with dissipation range. This 
was considered in Teufel & Schlickeiser 2002 and 2003 for 
pure slab geometry and in Shalchi & Schlickeiser 2003 for 
pure 2D and composite slab/2D geometry. 

In the current paper we compared the different contri- 
butions to the Fokker-Planck coefficient. We find that for 
a/i <C e the large scale contribution controlls the Fokker- 
Planck coefficient whereas for e <C a/i the small scale contri- 
bution is dominant. Thus for slab model turbulence quasilin- 
ear approximation is a good approximation except for CRs 
with pitch angle close to 90°. In this regime (fi m 0) non- 
linear effects can no longer be neglected (see Shalchi 2005). 
Because of the results of the current paper we come to the 
conclusion that previous calculations done in the quasilin- 
ear limit and for slab geometry are valid. In non-slab mod- 
els (e.g. the composite slab/2D model of Bieber et al. 1994) 
however, nonlinear effects always play a significant role, be- 
cause perpendicular diffusion itself can have a strong influ- 
ence on pitch-angle scattering (see Shalchi et al. 2004). 
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